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Abstract. We develop a systematic procedure of finding integrable "relativistic" (regular one- 
parameter) deformations for integrable lattice systems. Our procedure is based on the integrable 
time discretizations and consists of three steps. First, for a given system one finds a local dis- 
cretization living in the same hierarchy. Second, one considers this discretization as a particular 
Cauchy problem for a certain 2-dimensional lattice equation, and then looks for another mean- 
ingful Cauchy problems, which can be, in turn, interpreted as new discrete time systems. Third, 
one has to identify integrable hierarchies to which these new discrete time systems belong. These 
novel hierarchies are called then "relativistic", the small time step h playing the role of inverse 
speed of light. We apply this procedure to the Toda lattice (and recover the well-known relativistic 
Toda lattice), as well as to the Volterra lattice and a certain Bogoyavlensky lattice, for which the 
"relativistic" deformations were not known previously. 
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1 Introduction 

The theory of integrable differential-difference, or lattice, systems is by now a well developed 
and well understood subject. Nevertheless, some intriguing questions remain open, and the 
aim of this paper is to close one of them: what is relativistic Volterra lattice? Let us start 
with the necessary background. 

Certainly, two most celebrated and well studied integrable lattice systems are the Toda 
lattice (TL), 

b k = a k -a k -i, a k = a k (b k+1 - b k ) (1.1) 
and the Volterra lattice (VL), 

Uk = u k (v k - v k -i), v k = v k {u k+ i - u k ) (1.2) 

For readers more familiar with another forms of these systems, we recall that the Newtonian 
form of the Toda lattice, 

x k = e Xk + 1 ~ Xk - e Xk ~ Xk - 1 (1.3) 
equivalent also to the Hamiltonian one, 

x k = p k , p k = e Xk + 1 ~ Xk — e Xk ~ Xk - 1 (1.4) 

is recovered from ( |1 . 1|) , if the variables a k , b k are parametrized according to Manakov- 
Flaschka formulas 

b k = Pk} a k = e Xk +±- Xk (1.5) 
while the more convenient form of the Volterra lattice, 

a k = a k (a k+1 - a fc _ x ) (1.6) 

arises from (|1.2| ) upon the re-naming 

u k = a 2k -i, v k = a 2k (1.7) 

These two lattice systems are connected to one another in two different ways. On the one 
hand, the Volterra lattice ( |1.6|) is a restriction to the manifold b k = of the second flow of 
the Toda hierarchy. On the other hand (and this connection will be of a primary interest for 
us here) the flows (|1 . 1| ) and ( |1.2j ) are connected by a Miura map, or, better, by two different 
Miura maps M.i$ '■ (u,v) i— >■ (a, b): 

M (b k = u k + v k ^ M 2 :l bk = Uk + Vk (18) 
a k = u k v k ' \ a k = u k+ iv k 
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Some time ago a remarkable discovery was made by Ruijsenaars [[Rj] in the area of inte- 
grable lattice systems: he found a relativistic generalization of TL (RTL). The corresponding 
system may be viewed as a regular deformation of TL: 



d k = (1 + hd k )(a k - a k - 



dk — o,k(dk+i — dk + ha k+ i — ha k - 



1.9) 



(the reason for choosing here d instead of b will become clear in the main text). Under the 
parametrization 



dk = (ehpk _ i) / h) ak = e x k+ i -x k + hp k 
(which is a regular deformation of (|1.5|) ) the equations of motion ([I 

p k = e x k+i—x k + hp k _ e x k — Xk-i + hp k -i 

l + hx k = e h P k (l + h 2 e Xk + 1 ~ Xk 
which implies also Newtonian equations of motion 



(1.10) 

may be presented as 
(1.11) 



x k — (1 + hxk){l + hx 



[1 + /i.i fc _i)(l + hx k )- 



jXk %k—l 



1.121 



1 _|_ /^gXfc+i — x k v /K ' i _|_ h 2 e Xk ~ Xk ~ 1 

Mathematical structures related to RTL, including Lax representations, multi-Hamiltonian 
structure, and so on, were further investigated in |[BR|| , ||ZTOF|] , |Sl| , |S2| . A general ap- 
proach to constructing "relativistic" generalizations of integrable lattice systems, applicable 
to the whole lattice KP hierarchy, was proposed in ||GK1 ] . 



Paradoxically, up to now nobody seems to know what relativistic Volterra lattice is. We 
propose here an answer to this question. Let us stress that we are not concerned here with 
relativistic invariance. Instead, our aim is to construct integrable lattice systems, which 
are regular one-parameter deformations of VL and are Miura-related to RTL in the same 
manner as VL is related to TL. Needless to say that the "relativistic Miura maps" we are 
looking for have to be regular deformations of the standard ones. 

Surprisingly, these Miura maps can be choosen to be identical with the nonrelativistic 
ones. There exist two different systems that are sent to the RTL by either of two Miura 
maps ( p..8|) , .Mi^'M, v) = (a, d), so that "the relativistic Volterra lattice" actually splits into 
two systems: 

ii k = u k (v k - v k _ x + hu k v k - huk-iv^) 

(1.13) 

Vk = v k {u k+l -u k + hu k+1 v k+1 - hu k v k ) 



and 



u k (v k - Ufc-i + hu k+ iv k - hu k v k -i) 
v k {u k+ i -u k + huk+iVk - hu k v k -i, 



;i.i4) 
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We suspect that our way to finding these systems might be more significant than the sys- 
tems themselves. Namely, our route is through the theory of integrable time discretizations. 
It was discovered in |ST[ that certain integrable discretization of the usual TL, 

(1.15) 



b k = b k + h(a k - o fc _i), o fc (l + hb k ) = a k (l + hb k+1 ) 



shares the integrals of motion and the invariant Poisson structure with RTL (upon the 
change of variables b k = d k + ha k _i). In other words, this discretization belongs to the RTL 
hierarchy. In [ |PGK|| , ||S3|| it was noticed that this discretization is connected to other one, 
belonging to the TL hierarchy: 



b fc = b fc + h(& k - a fe _ 



a fe (l + hh k ) = a fc (l + hb k+1 ) 



;i.i6) 



These two 1+1-dimensional discrete systems (with one discrete space coordinate and one 
discrete time) may be seen as resuling from one and the same 2-dimensional discrete system 
(living on a 2-dimensional lattice) by posing an initial value problem in two different ways. 
In other words, the Cauchy data for two discretizations are prescribed on two different 
"discrete curves" on the 2-dimensional lattice. Being very close on the 2-dimensional lattice, 
the equations (|1.15|) , ( |1.16|) have nevertheless very different properties, such as invariant 
Poisson brackets, integrals of motion, Lax matrices, etc. - all in one, they belong to different 
hierarchies. 

Thus, our strategy to finding the "relativistic" deformations for continuous time lattice 
systems (with one discrete space coordinate) may be described as follows. 

First, for a given integrable lattice system, one constructs integrable time discretization 
belonging to the same hierarchy. Such time discretizations appeared first in |[AL|| , [ GK2 1 , 
and a systematic procedure was developed in fS"3| , @. As a rule, this approach results 



in nonlocal equations of motion. However, these nonlocal discrete time equations of motion 
may be brought into a local form with the help of the so called localizing changes of variables, 
which were found in ||S6|| for a large set of examples. 

Second, and this is a crucial step in finding new hierarchies, one considers the result- 
ing discrete time system as a system on a two-dimensional lattice, and tries to find new 
meaningful initial value problems for this two-dimensional lattice. This is close in spirit to 
the constructions in |[PNC|| , |[NPCQ|| , [PN|| . The resulting discrete time systems belong to 
integrable hierachies distinct from the original one (this fact is not stressed in the papers 
just mentioned). 

As the third and final step, one has to identify these novel hierarchies, in particular, to 
find integrals of motion, invariant Poisson structures, and higher flows. 

In this paper we first recall how this program could be realized to find RTL (although 
the actual way to discovering this system was quite different), and then demonstrate how 
the relativistic Volterra lattice may be derived. In particular, we show that ( |1.13| ) and ( |1.14| ) 
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are the simplest flows of the hierarchies, to which the following two explicit discretizations 
of VL belong: 

u k (l + hv k -i) = u k (l + hv k ), v k (l + hu k ) = v k {\ + hu k+1 ) (1.17) 

and 

u k (l + fovfc-i) = u fc (l + hv k ), v k (l + hu k ) = v k (l + hu k+x ) (1.18) 

respectively. One can see that the corresponding constructions may be described as a fac- 
torization of RTL, in a complete analogy with the nonrelativistic case. Finally, we show how 
to generalize these results to some of the Bogoyavlensky lattices. 



2 General framework 

2.1 Lax equations and representations 

Our approach to integrable lattice systems is based on the notion of Lax representations. 
We consider Lax equations of one of the following types: 



L,7T + (f{L)) 



(2.1) 



or 



U))-Mf(T 3 )).L 



Lj • 7T_ 



+ 7T_ 



(mi) ■ Lj (2.2) 



Let us explain the notations. 

Let g be an associative algebra. One can introduce in g the structure of Lie algebra in a 
standard way. Let g + , g be two subalgebras such that as a vector space g is a direct sum 
g = g + ® g_. Denote by 7r ± : g \— > g ± the corresponding projections. Finally, let / : g i— > g 
be an Ad-covariant function on g, and let L stand for a generic element of g. Then ( [2.1|) is 
a certain differential equation on g. 

Further, let g = <S>JLi 9 be a direct product of m copies of the algebra g. A generic 
element of g is denoted by L = (L 1; . . . , L m ). We use also the notation 

Tj — Tj\h) — Lj • . . . ■ L\ • L m ■ . . . ■ Lj + i (2-3) 



Then ( |2.2| ) is a certain differential equation on g. Such equations are sometimes called Lax 
triads. 

One says that (|2.1|) , resp. ( [2.2|) , is a Lax representation of a Hamiltonian flow 



x = {H, x} 



(2.4) 



4 



on a Poisson manifold \X, {•, -}J, if there exists a map L : X i— > gt (resp. L : A? i— > g) such 
that the former equations of motion are equivalent to the latter ones. Let us stress that 
when considering equations ( |2.1|) , resp. (|2.2|) in the role of Lax representation, the letter 
L (resp. L) does not stand for a generic element of the corresponding algebra any more; 
rather, it represents the elements of the images of the maps L : X \— > g and L : X i— > g, 
correspondingly. The elements L(x), resp. L(x) (and the map L, resp. L, itself) are called 
Lax matrices. 

2.2 r— matrix Poisson brackets 

Recall that there exist several constructions of Poisson brackets on associative algebras im- 
plying the Lax form of Hamiltonian equations of motion. We recall here some of them. 

Suppose that g carries a nondegenerate scalar product (•,-), bi-invariant with respect to 
the multiplication in g. Let R be a linear operator on g. 

Definition 2.1 ||STS[ | A linear r-matrix bracket on g corresponding to the operator R is 
defined by: 

{<p, ^h(L) = I( [R(V^(L)), + [Vv?(L), R(W(L))] , L ) (2.5) 

// this is indeed a Poisson bracket, it will denoted by PB^R). 



Theorem 2.2 [5T5 A sufficient condition for (|2.5|) to define a Poisson bracket is given by 



the modified Yang-Baxter equation for the operator R, mYB(R; a), which reads 

[R(w),R(t>)] -R([R(m),w] + [u,R{v)\) = -a[u,v] \fu,veg (2.6) 
Now let Ai, A 2 , S be three linear operators on g, A± and A 2 being skew-symmetric: 

Al = -Ai, A; = -A 2 (2.7) 



Definition 2.3 [ 52 ] A quadratic r-matrix bracket on g corresponding to the triple A±, A 2; 
S is defined by: 

W^h(L) = ^(A x {d'^{L)),d^{L))-^(A 2 {d V {L)),d^{L)) 

+\ (S{ML)), d'^(L) ) - \ (S*(dV(i)), dif>(L) ) (2.8) 

where we denote for brevity 

dip(L) = L-V^(L), d'ip(L) =V<p(L) ■ L (2.9) 
If this expression indeed defines a Poisson bracket, we shall denote it by PB 2 (Ai, A 2 , S). 
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In what follows we shall usually suppose the following condition to be satisfied: 

Ai + S = A 2 + S* = R (2.10) 

Then a linearization of PB 2 (A 1; A 2 , S) in the unit element of g coincides with PB^R), and 
we call the former a quadratization of the latter. 



Theorem 2.4 | S2 | A sufficient condition for (2^) to be a Poisson bracket is given by the 



equations ( |2.10| ) and 

mYB(R; a), mYB(A i; a), mYB(A 2 ; a) (2.11) 

One of the most important properties of the r-matrix brackets is the following one. 

Theorem 2.5 Ad-invariant functions on g are in involution with respect to the bracket 
PBx(R) and with respect to its quadratizations PB 2 (A 1; A 2 , S). The Hamiltonian equations 
of motion on g corresponding to an Ad-invariant Hamilton function if, have the Lax form 



L = - 
2 



L,R(/(L)) (2.12) 



where f(L) = Vy(L) for the linear r-matrix bracket, and f(L) = df(L) for its quadratiza- 
tions. 

Quadratic r-matrix brackets have interesting and important features when considered 
on a "big" algebra g = 0^=19- This algebra carries a (nondegenerate, bi-invariant) scalar 
product 

m 

((L,M» = £(L fc ,M fc > 

k=l 

Working with linear operators on g, we use the following natural notations. Let A : g i— > g 
be a linear operator, let ^A(L)) be the ith component of A(L); then we set 



A(L)) = £( A k(^') (2.13) 

3=1 

For a smooth function $(L) on g we also denote by Vj$, c/j$, c/j$ the j'th components of 
the corresponding objects. 

Now let Ai, A 2 , S be linear operators on g satisfying conditions analogous to ( |2.7| ) and 
to ( [2. 11| ). One has, obviously: 



A^)* = -(Ai)ji, ((A 3 )y)* = -(A a )i<, (Sy)' = (S*) i4 
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Then one can define the bracket PB 2 (A 1; A 2 , S) on g. In components it reads: 

-r m i m 

{<&,*} 2 (u) = - E {(AMdpwp)-- E {{AMd&tdi* 



i m 

+2 E (SyCd^.d/* 



i m 

E ((s*) y (d;$),4* 



2 .4 , 



(2.14) 



Theorem 2.6 ||S5|| Let g be equipped with the Poisson bracket PB 2 (A 1; A 2 , S). Suppose that 
the following relations hold: 



(Ai) i+ ij + i 



H+i,j 



(A 



2)i,j 



for i ^ j 



( A i)i+i,i+i - + (S) i+ i,j - (S*) iii+ i = for all l<j<m 



Then each map Tj : g i— > <7 (|2.3| ) zs Poisson, if the target space g is equipped with the Poisson 
bracket 

PB 2 ((Ai) i+ i iJ+1 , (A 2 )jj, (S)j + ijJ 

Hamilton function of the form $(L) = tp(L m ■ ... ■ L\), where ip is an Ad-invariant function 
on g, generates Hamiltonian equations of motion on g having the Lax form: 



BjLj, 



-R,(^(T,)) 



(2.15) 



where 



R, 



(Ax! 



+ (S) 



3+1,3 



A 2 )« + (S*) 



(7n a// the formulas the subscripts should be taken (mod m)). 



We have discussed above the r-matrix origin of Lax equations. If one is concerned with 
a Lax representation of a Hamiltonian flow (|2.4j) on a Poisson manifold {X, {•,-}), then 
finding an r-matrix interpretation for it consists of finding an r-matrix bracket on g (or on 
g) such that the Lax matrix map L : X i— > g (resp. L : ^ i— > g) is a Poisson map. 



2.3 Factorizations and integrable discretization 



A further remarkable feature of the equations ( |2.1| ) and fl2.2| ) is a possibility to solve them 
explicitly in terms of a certain factorization problem in the Lie group G corresponding to g 
|Syfl , [ |5TS|| , [|RSTS|| . (Actually, this can be done even in a more general situation of hierarchies 
governed by .R-operators satisfying the modified Yang-Baxter equation, see |[RSTS|| ). The 
factorization problem is described by the equation 



u = u + {u) -n_(to, ueg, n±(u)eG± 



(2.16) 
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where G± are two subgroups of G with the Lie algebras g±, respectively. This problem has 
a unique solution in a certain neighbourhood of the group unit. 

In the situations we consider in the sequel G will be a matrix group, and we write the 
adjoint action of the group elements on g as a conjugation by the corresponding matrices. 
In this context we write Il^ 1 (?7) for (II±({7)) _1 . Correspondingly, we call Ad-covariant 
functions g i— > g also "conjugation covariant". This notation has an additional advantage of 
being applicable also to functions g i— > G. 

Based on the above mentioned explicit solution, the following recipe for integrable dis- 
cretization was formulated in In all the difference equations below we use the tilde 
to denote the discrete time shift, so that, for example, L = L(t + h), if L = L(t), t e hZ. 

Recipe. Suppose you are looking for an integrable discretization of an integrable system 
fl2.4| ) allowing a Lax representation of the form ( |2.1| ) . Then as a solution of your task you 
may take the difference equation 

L = n+\F(L)) ■ L ■ U+(F(L)) = n_(F(L)) • L ■ UZ^L)) (2.17) 

with the same Lax matrix L and some conjugation covariant function F : g i— > G such that 

F(L)=I + hf(L)+0(h 2 ) 



Analogously, if your system has a Lax representation of the form ( |2.2j ) on the algebra g ; then 
you may take as its integrable discretization the difference Lax equation 

L j = ITT^FCr,)) • Lj • n + (F(Tj_i)) = II_(F(T,)) • L 3 • ir 1 ^^)) (2.18) 

with F as above. 

Obviously, by construction, the discretizations obtained in this way share the Lax ma- 
trix and therefore the integrals of motion with their underlying continuous time systems. 
Moreover, they share also the invariant Poisson brackets. Indeed, the above mentioned 
factorization theorems imply that the maps ( p. 17 ), Q2.18 ) are the time h shifts along the 



trajectories of the corresponding flows ( |2.1| ), ( |2.2| ) with 

f(L) = h- 1 \og(F(L)) = f(L) + 0(h). 

This, in turn, implies that if all flows of the hierarchy ( |2.1| ) [resp. ( |2.2| )] are Hamiltonian with 
respect to a certain Poisson bracket, then our discretizations are Poisson maps with respect 
to this bracket. In particular, if the Lax matrices L [resp. L] form a Poisson submanifold for 
some r-matrix bracket, then this submanifold is left invariant by the corresponding Poisson 
map ( |2.17| ) [resp. ( |2.18 )]. We shall say that our recipe yields discretizations living in the 



same hierarchies as the underlying continuous time systems. 
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2.4 Basic algebras and operators 

In what follows we fix an algebra g suited to describe various lattice systems with periodic 
boundary conditions, as well as certain operators taking part in the corresponding r-matrix 
constructions. So, starting from this point the symbols g, g±, G, G±, tt±, U±, R, Ai, A 2 , S 
will carry only the following meanings. 

The algebra g is a twisted loop algebra over gl(N), consisting of Laurent polynomials L(\) 
with coefficients from gl(N) and a natural commutator [u\^,v\ k ] = [u,v]\ : > +k , satisfying an 
additional condition 

fiL(A)^ 1 = L(uX), where Q = diag(l, . . . , c/^ 1 ), u = exp(2m/N) 

In other words, elements of g have the following structure: 

L(\)=Y J Z {P) £ P (2-19) 

v 

In this formula = diag (ii , . . . , 1% ) are diagonal matrices, and 

N N 

£ = A Ek+i^k, £ = A ^2 E^k+i 

k=l k=l 

are the shift matrices (here and below Ej^ stands for the matrix whose only nonzero entry is 
on the intersection of the jth row and the kth column and is equal to 1; we set £/jy+i,Ar — -^i,AT) 
En,n+i = ^at,1) and in general understand all subscripts (mod N)). 
The nondegenerate bi-invariant scalar product on g is chosen as 

(L(A), M(A)) = tr(L(A) • M(A)) (2.20) 

the subscript denoting the free term of the formal Laurent series. This scalar product 
allows to identify g* with g. 

The two subalgebras g + , g_ such that g = g + ® g are choosen as 

g+= lj2i {p) £ p \, g- = \ y Z i(p)£P \ ( 2 - 21 ) 

[p>0 J [p<0 J 

The corresponding decomposition L = 7t+(L) + 7r_(L) of an arbitrary element L E g will be 
called the generalized LU decomposition. 

We denote also by g the commutative subalgebra of g consisting of A-independent 
diagonal matrices (i.e. of matrices £^£°). The linear operator on g assigning to each 
element L(X) ( |2.19| ) its free term £^ will be denoted by Pq. 
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The group G corresponding to the twisted loop algebra g i is a twisted loop group, con- 
sisting of GL(iV)-valued functions U(X) of the complex parameter A, regular in CP 1 \{0, 00} 
and satisfying QU^X)^ 1 = U(uX). Its subgroups G + and G_ corresponding to the Lie 
algebras g + and g_, are singled out by the following conditions: the elements U(X) E G + 
are regular in the neighbourhood of A = 0, and the elements U(\) G G_ are regular in 
the neighbourhood of A = 00 and take in this point the value U(oo) = I. We call the 
corresponding n + n_ factorization the generalized LU factorization. 

The basic operator governing the hierarchies of Lax equations, is: 

R = 7T_|_ — 7T_ (2.22) 

Denote by Ro, Po its skew-symmetric and its symmetric parts, respectively: 

Ro = (R-R*)/2, P = (R + R*)/2 

Obviously, this definition of Po is consistent with the previous one. 
Let the skew-symmetric operator W act on g according to 

k<j k>j 

and on the rest of g according to W = W o P . Finally, define: 

A 1 = R + W, A 2 = R -W, S = P -W, S* = P + W (2.23) 

These operators will be basic building blocks in the quadratic r-matrix brackets appearing 
below. 

3 Reminding the Toda lattice case 
3.1 TL 

We consider the equations of motion of TL fll.ip under periodic boundary conditions: all 
subscripts are taken (mod N). The phase space of this system is 

T = M 2N {b l ,a 1 ,...,b N ,a N ) (3.1) 

(recall that a = ajv, b N+1 =61). 

There exist three compatible local Poisson brackets on T such that the system TL is 
Hamiltonian with respect to each one of them, see |[K1|| . We adopt once and forever the 
following conventions: the Poisson brackets will be defined by writing down all nonvanishing 
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brackets between the coordinate functions; the indices in the corresponding formulas are 
taken (mod N). 

The "linear" Poisson structure on T is defined by the brackets 

\bk, Qfc}i = — dk, {a k , &fc+i}i = — a>k (3.2) 
the corresponding Hamilton function for the flow TL is given by: 

j JV N 

H 2 (a,6) = -E^ + E^ (3-3) 

z k=l k=l 

The "quadratic" Poisson structure has the following definition: 

\bk,Q'k}2 = — a kbk, { a k,bk+i}2 = ~ a kbk+i 

(3.4) 

{bk,bk+l}2 = - Qfcj { a k, a k+l}2 = — a k+l a k 

The Hamilton function generating TL in this bracket is: 

N 

R 1 (a,b) = Y,bk (3.5) 
fc=i 

Finally, the "cubic" bracket on T is given by the relations 

{bk,a k }z = -a k (b 2 k + a k ), {a k ,b k+1 } 3 = -a k {b\ +1 + a k ), 

{bk,bk+i}3 = — &k{bk + bk+i), {ak,ak+i}3 = —2a k a k+ ib k+ i, (3.6) 

{bk,0-k+l}3 — — a k a k+l, {a>k,bk+2}3 — —O-kO-k+l 

The corresponding Hamilton function of the flow TL is: 

1 N 

HoM) = -$>gK) (3.7) 
A k=i 



The Lax representation of the Toda lattice JF), [M] lives in the algebra g introduced in 



the previous section. We shall work with the Lax matrix 

N N N 

L(a, b, A) = A E k+i,k + Y b k E kyk + A -1 Y a-kE k)k+ i = £ + b + aS' 1 (3. 

k=l k=l k=l 

where diagonal matrices a = diag(ai, . . . , ajv), b = diag(&i, . . . , b N ) are introduced. 
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The equations of motion ( |1 . 1|) are equivalent to the Lax equations 



L=[L,B]=-[L,A] (3.9) 

with 



N N 



B = tt + (L) = Xj2E k+ i, k + J2 b kE k , k = S + b (3.10) 

k=l k=l 
N 

A = vr_(L) = A -1 ^2 a k E k>k+ i = aS- 1 (3.11) 



k=l 



where 7r ± : g i— > g ± are the projections to the subalgebras g ± . 

Spectral invariants of the Lax matrix L(a, b, A) serve as integrals of motion of this sys- 
tem. Note that all Hamilton functions in different Hamiltonian formulations belong to these 
spectral invariants. For instance, 

H a (o,6) = i(trL 2 (a,6,A)) Q , E 1 (a,b) = (tr L(a, b, A)) q 

where the subscript " 0" is used to denote the free term of the corresponding Laurent series. 

All spectral invariants turn out to be in involution with respect to each of the Poisson 
brackets ( |3.2|) , (|3.4f) , ( |3.6|) . Most directly it follows from the r-matrix interpretation of the 
Lax equation Q3.9p. 

Theorem 3.1 The Lax matrix map L(a, b, A) : T i— > g is Poisson, if T carries {•, -}i and g 
carries PB!(R), and also if T carries {•, -}2 and g carries PB 2 (A 1; A 2 , S). 

The first statement is from [ AM ], the second one - from [|52"| . For other versions of such 
statements (including the r-matrix interpretation of the cubic bracket) see |DLT1|| , [pR]| , 
[MP]. 



3.2 TL -> dTL 

In order to find an integrable time discretization for the flow TL, we apply the recipe of the 
previous section with F(L) = I + hL, i.e. we take as a solution of this problem the map 
described by the discrete time Lax equation 

L = B l LB = ALA' 1 with B — n + (J + hL), A = n_(J + hL) (3.12) 
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Theorem 3.2 ||S3|| (see also ||GK2 |). Consider the change of variables T(a, b) i— > T(a,b) 
defined by the formulas 



h = h k + ha, 



■fc-i) 



a k = a fe (l + hh k ) 



(3.13) 



The discrete time Lax equation ( |3. 12| ) zs equivalent to the map (a, b) i— > (a, 6) which in the 
variables (a, b) is described by the following equations: 



bfc = b fc + /i(a fc - a fc _ 



a*(l + /ibfc) = a*(l + Kb, 



k+l, 



(3.14) 



Proof. The tri-diagonal structure of the matrix L implies the following bi-diagonal structure 
for the factors B, A: 



N N 

B = !!+(/ + hL) = 5^(1 + hh k )E Kk + h\Y J E l 



k+l,k 



k=l 



k=l 



N 



k,k+l 



(3.15) 



(3.16) 



fc=i 



The formulas ( |3.13| ) represent the matrix equation / + hT = BA which serves as a def- 
inition of the matrices B, A. Obviously, these formulas define a local diffeomorphism 
T(a, b) i — > T(a,b). Now, the Lax equation fl3.12|) is easy to see to be equivalent to 
I + hf = B 1 (BA)B = A(BA)A" 



or 



B A = A B 



(3.17) 



The equations of motion ([3.14 ) are now nothing but the componentwise form of the latter 
matrix equation. I 



The map (|3.14j) will be denoted dTL. Obviously, it serves as a difference approximation 
to the Toda flow TL (|1.1|) . The construction assures numerous positive properties of this 
discretization: in the coordinates (a, b) the map dTL is Poisson with respect to each one 
of the Poisson brackets ( |3.2| ), (|3.4j) , ( |3.6| ), it has the same integrals of motion as the flow 
TL, etc. Going to the coordinates (a, b) deforms the integrals of motion and the Poisson 



brackets. Moreover, since the inverse to the map ( |3. 13| ) is nonlocal, the invariant Poisson 
brackets become, generally speaking, also nonlocal in the coordinates (a, b). Remarkably, 
there turn out to exist such linear combinations of the basic invariant Poisson brackets whose 
pull-backs to the coordinates (a, b) are described by local formulas. 



Theorem 3.3 a) The pull-back of the bracket 

{vh + Mv}a 



(3.1i 
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on T(a,b) under the change of variables (|3.13| ) is the following bracket on T(a, b): 

{b fc ,a fe } = -a*(l + hh k ), {a k , b k+1 } = -a fc (l + hb k+1 ) (3.19) 

b) The pull-back of the bracket 

{v} 2 + Mv}a (3-20) 
on T(a,b) under the change of variables (|3.13| ) is the following bracket on T(a, b); 

{b fc , a fc } = -a k (b k + ha k )(l + hb k ), {a k , b k+1 } = -a fc (b fc+ i + ha k )(l + hb k+1 ) 
{b k , b k+1 } = -a k (l + hb k )(l + hb k+1 ), {a k , a k+1 } = -a fc a fc+1 (l + hb k+l ) 

(3.21) 

c) The brackets ( |3.19| ), G3.21| ) are compatible. The map dTL ( |3.14| ) is Poisson with respect 
to both of them. 

Proof. To prove the theorem, one has, for example, in the (less laborious) case a) to verify 
the following statement: the formulas (|3.19| ) imply that the nonvanishing pairwise Poisson 
brackets of the functions ( |3.13| ) are 

{b k , a k } = -a k (l + hb k ), {a k , b k+1 } = -a k (l + hb k+1 ) 

{b k ,b k+ i} = —ha k , {ak, a k+i} = ~ha k a k+ i 

This verification consists of straightforward calculations. I 

The map ( p,14[ ) was first found in [|HT1|| , along with the Lax representation, but without 
discussing its Poisson structure and its place in the continuous time Toda hierarchy. 

We can now determine the hierarchy of continuous time lattice equations to which the 
map ( |3.14| ) belongs. Clearly, this is the Toda hierarchy pulled-back under the map ( |3.13| ). 
The previous theorem allows to calculate the corresponding equations of motion in a sys- 
tematic (Hamiltonian) fashion. 

Theorem 3.4 |[K1|| . The pull-back of the flow TL under the change of variables ( |3.13| ) is 
described by the following differential equations: 

b fe = (a fc - a fc _i)(l + hb k ), a k = a fc (b fc+ i - b k ) (3.22) 

Proof. To determine the pull-back of the flow TL, we can use the Hamiltonian formalism. 
An opportunity to apply it is given by the Theorem |3.3| . We shall use the statement a) only. 
Consider the function /i _1 H 1 (a,6) = h^ 1 J2 k= ib k - It is a Casimir of the bracket {■, and 
generates exactly the flow TL in the bracket h{-, -}2. Hence it generates the flow TL also in 
the bracket (|3.18|) . The pull-back of this Hamilton function is equal to h^ 1 Y, k =i(bk + ha k ^i) . 
It remains only to calculate the flow generated by this function in the Poisson brackets ( |3.19| ). 
This results in the equations of motion Q3.22|) . ■ 
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3.3 dTL -> explicit dTL 



Now consider the equations (|3.14j) as equations on the 2-dimensional lattice. In other words, 
we attach the variables = bfc(t), = a&(i) to the lattice site (k,t) G Z x hz. A 
linear change of independent variables (k, t) i— ► r) = t + fc/i) mixes space and time 
coordinates, which is equivalent to changing the Cauchy path on the lattice to {r = 0} = 
{t = —kh}. For another instances of such staircase (or sawtooth) Cauchy paths on the 
lattice, see fPNO, [NPUO , [ |FV| |. To deal with the new initial-value problem, denote 

b k (r) = h k (t) =b k (r - kh), a k (r) = a fe (t) = a fc (r - kh) (3.23) 

Denoting the /i shift in r still by the tilde, it is easy to see that the variables a k = a k (r), 
b k = b k (r) satisfy the following difference equations: 



b k = b k + h(a k - a k _ 



a k (l + hb k ) = a k (l + hb k+1 ) 



(3.24) 



These equations serve as an explicit discretization of the flow TL. Indeed, they allow to 
calculate (a, b) explicitly, if (a, b) is known (first b, then a). 

Of course, it is tempting to conclude about the integrability of the system ( |3.24j ) from the 
known integrability of the equations (|3.14j ). However, we do not know any general statements 
allowing such conclusions, so that the integrability of (|3.24j) has to be proved independently. 
Moreover, we are not aware of any systematic procedure allowing to determine invariant 
Poisson brackets of the map ( |3.24 ) starting with the known invariant Poisson brackets of 
the map ( |3.14| ) . Hence even the definition of integrability is nontrivial in the case of ( |3.24| ) . 
Our approach to these problems will be based on Lax representations. 

Recall that the Lax equations may be considered as compatibility conditions of linear 
problems. In particular, the discrete time Lax equation ( |3.17| ) has two versions: 



BA = A(BA)A 



-i 



B 1 (BA)B 



The first one is a compatibility condition of two linear problems: 

BAip = fj,tp, $ = Aip 
while the second one is a compatibility condition of 

BAip = fit/;, B~V 



(3.25) 



(3.26) 



Here ip = (ipi, . . . , is an auxiliary vector, fi is a spectral parameter (having nothing 
in common with the spectral parameter A entering the elements of the algebra g such as 
L(a,b,X)). 
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These two pairs of linear problems are equivalent, correspondingly, to the following ones: 



Btp = fix/;, A0 = -0 (3.27) 

and 

Aip = B^p = ip (3.28) 

We are now in a position to derive from these two pairs of auxiliary linear problems two 
different Lax representations for the map Q3.24|) . To this end, we assume that the variables 
ip k (t) are also attached to lattice sites (k,t), and perform in the above equations the change 
of variables (k,t) i— > (k,r) = (k,t + kh). We denote, in addition to the variables a k , b k 
introduced above, also 

4>k( r ) =Mt) =ipk(r-kh) 



Theorem 3.5 The map ( 3.24 ) allows the following Lax representation: 

L = D » X L D = C l L C DC 1 = C l D (3.29) 

with the Lax matrix 



L = 

where 



I N N N \ 

DC~ 1 -I)/h =[XJ2 E k +i,k + E( b fc - hak-jEuc + A" 1 £ a k E Kk+1 CT 1 (3.30) 

\ k=l k=l k=l ) 



N N 



D = J2( 1 + hb k-h 2 a k ^ 1 )E kk + h\J2E k+1 , k (3.31) 

k=l k=l 
N 

C = I-h\- l Y,a k E k>k+1 (3.32) 



k=l 



Proof. The statement may be, of course, easily verified, but we want to show how it can be 
derived. To this end we start with ( |3.27[ ). In coordinates these equations read: 

(1 + hb k )ipk + hX^k-i = (iif) k , ^pk = ^k + hX~ 1 a. k ip k+1 

This implies, after the change of variables (k,t) i— > (k,r): 

(1 + hb k )$ k + hXij)^ = iMil> k , i> k = il> k + hX^a^^ 

After simple manipulations the latter system may be put into the form 

(1 + hb k - /i 2 a fc -i)V> fe + hXtl)^ = ni> k , ip k = ^ k - /jA -1 a fc $ fc+1 
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In matrix notations the latter equations look like 

with the matrices D, C introduced above. The compatibility condition of these two linear 



problems coincides with ( |3.29|) . 



Theorem 3.6 The map ( |3.24|) allows the following Lax representation: 

L + = A- 1 J r -L + -J 7 ' 1 A 

with the matrices 



(3.33) 



/ N N N \ 

L+ = T~ x ( A ^2 E k+i,k + E ( b k - ha k ^)E kk + A" 1 E a k E k>k+1 (3.34) 

V k=l k=l k=l I 



N 



f = I - h\J2 Ek+i, k = I - hS 

k=l 

A = diag(l + hbx, . . . , 1 + hb N ) 



(3.35) 
(3.36) 



Proof. We perform with ( |3.28| ) transformations analogous to those of the previous proof. 
In coordinates (|3.28| ) reads: 

ip k + /iA _1 a fc ^fc+i = H4>k, (1 + hb k )tp k + hAV> fe _i = ip k 

This implies, after the change of variables (k,t) i— > (k,r): 

if) k + hX^a^^ = /i^ fc , (1 + hb k )^ k + hXifj^ = <0 fc 

Straightforward manipulations imply: 

( 1 + hb k - h 2 a k 1 + ^ fc ) ip k + hX^ak 1 + ^ k ip k+1 = fi(^ k - hXijj^), 
\ 1 + hb k+1 J 1 + hb k+1 

1 



1 + hby 



This may be put in a matrix form: 
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with the matrices J 7 , A as above, and 



P=^ l + : ■ r r — E kk + h\ }_^a k — E k , k+1 (3.37) 

k=i V 1 + "■"fc+i / fc=i 1 + ftt, fe+i 

Now the compatibility condition of the latter two linear problems reads: 

T~ X P = A' 1 J 7 ■ T~ X P ■ J 7 ' 1 A 



In order to bring this into the required form (|3.33 ), notice that the equations of motion 
( |3.24| ) imply the following simple formula: 

N N 
P = J2{1 + hb k - h 2 a k _ x )E kk + h\- 1 a k E k>k+1 (3.38) 

k=l k=l 

so that 

jr-ip = I + hL+ 

This completes the proof. ■ 

3.4 Explicit dTL RTL 



Now that we have two different Lax representations for the map (|3.24 ), with the Lax matrices 



L±, the natural question is the one about their relation. It is easy to see that this relation 
may be described as follows: 

where a = (a 1 a 2 ...a N ) 1 / N and 

Oi = diag(l, a _1 ai, a~ 2 aia 2 , . . . , a _Ar+1 aia 2 ...aAr_i) 



Both Lax matrices for the explicit dTL (|3.24 ) are substantially different from the one for the 



implicit dTL ( |3.14|) (the latter Lax matrix is the usual Toda one (|3.8| ) with the variables (a, b) 
parametrized according to ( 3.13Q ). The next natural question is: what is the continuous time 



hierarchy, to which ( |3.24 ) belongs? In other words, how can one get continuous time flows 
with the Lax matrix L±, and what is the underlying invariant Poisson structure(s)? 
The answer is known for both Lax formulations. Introduce the phase space 

n = R 2N {d l ,a l ,...,d N ,a N ) (3.39) 
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Theorem 3.7 ||S2|| The Lax matrix map L_(a, d, A) : 1Z i— > g, where 



/ N N N \ 

L = (DC 1 - J)//i = A ]T £ fc+ i, fe + £ d fe £ fcfc + A" 1 £ a fc £ fcifc+1 C 1 

\ k=l k=l k=l ) 



(3.40) 



N N 

D(a,d,X) = ^2(l + hd k )E kk + hXj2E k +i,k 

k=l k=l 



N 



C(a, d, A) = I -h\ 1 J2 a k E k)k+ i 

k=l 



(3.41) 
(3.42) 



is Poisson, if TZ is equipped with the Poisson bracket 

{d k , a k }i = — a k , {a k ,d k+ i}i = —a k , {d k , d k+ {\i = ha k (3.43) 
and g is equipped with PBi(R), and also iflZ is equipped with the Poisson bracket 
{d k ,a k } 2 = —a k d k , {a k ,d k+ i} 2 = —a k d k+ \ 

{d k ,d k+ i} 2 = —a k , {a k ,a k+ i} 2 = —a k a k+1 



(3.44) 



and g is equipped with PB 2 (A 1; A 2 , S) . The hierarchy of continuous time flows is of the usual 
form 



±7T±(/(L_)) 



(3.45) 



The "first" flow of the hierarchy, corresponding to f(L) = L, coincides with RTL ( |1.9|) and 
allows a Lax representation with the auxiliary matrix 



N N 

n+{LJ) = ^2(d k + ha k ^i)E kk + A ^ E k+1)k 

k=l k=l 



(3.46) 



The map 



L = D X L D = C l L_C 

is interpolated by the flow of this hierarchy with f(L) = h^ 1 log(J + hL). 



Actually, an additional information is available from fj'2\ . The evolution of the factors C, 
D for the flows of the RTL hierarchy (|3.45 ) is described by the Lax triads 



D = ±D Ti ± {f{C- l D)) T MfiDC- 1 )) ■ D 
C = ±C- n ± {f{C- l D)) T n ± (f(DC- 1 )) ■ C 
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In particular, for the flow RTL (|1.S| ) the corresponding auxiliary matrices in the Lax triads 
are ( |3.46p and 



N N 

v + ((C- l D - I)/h) = J2(dk + ha k )E kk + A £ E k+1 , 



(3.47) 



k=l 



k=l 



For the case of the quadratic bracket an r-matrix interpretation in g®g is possible. Namely, 
the Lax matrix map (£>, C _1 ) > g = g®g is Poisson, if 7Z is equipped with the bracket 
{v}i + Mv>2, i-e. 

{d k ,a k } = -(l + hd k )a k , {a k , d k+1 } = -(1 + hd k+1 )a k , {a k ,a k+1 } = -ha k a k+l 

(3.48) 

and g ® g is equipped with /iPB 2 (Ai, A 2 , S) with the operators 



Ax 



Ai -S 
S* Ai 



A 9 



A 2 -S* 
S A 2 



S S 
S -S* 



(3.49) 



Theorem |2.6| then assures that the monodromy matrix maps DC and C X D are Poisson, 
if the target g is equipped with /iPB 2 (Ai, A 2 , S), which implies the Poisson property for the 
Lax matrix maps L = (DC^ 1 — I)/h and (C^D — I)/h, if the target g is equipped with 
PB 1 (R) + /iPB 2 (A 1 ,A 2 ,S). 

For the hierarchy attached to the Lax matrix L + , somewhat less detailed information is 
available. 



Theorem 3.8 |BR|| , ||GK1|| The hierarchy of Lax equations related to the Lax matrix 

/ N N N \ 

L+(a, d, A) = (J - hEy 1 A ^ E k+1)k + ^ d k E kk + A" 1 ^ a k E k>k+l 

\ k=l k=l k=l ) 



(3.50) 



is given by the formula 



L + ,±7r ± {f{L + ))+a(f(L + )) 



where a is a linear map from g into the set of diagonal matrices defined as 

P p<0 

In particular, the "first " flow of this hierarchy corresponding to f{L) = L coincides with 
RTL (|1.9|) and allows the Lax representation 



A, L A 
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with the auxiliary matrix 



N N 

A = n-{L + ) - a{L + ) = A -1 a k E k:k+1 - h ^ a k E kk 

k=l k=l 

It has to be mentioned that each of the formulations above has its specific advantages. 



Theorem |3.7| has an advantage of including the RTL hierarchy into the standard framework 
of the lattice KP hierarchy. On the other hand, Theorem |3.8|, found for RTL in ||BR||, has 



an advantage of being a particular case of a much more general construction, due to ||GK1 



which delivers a "relativistic" generalization of the whole lattice KP hierarchy. An r-matrix 



interpretation of the latter result is unknown (see, however, |PR|| for a linear bracket in the 
open-end case). 

So, upon the identification 

d k = b k — ha k -i 

the explicit discrete time Toda map ( p. 24 ) belongs to the relativistic Toda hierarchy. 
Remarks. 

1. In the coordinates (a, b) the equations of motion Q1.9Q take the form 

b k = a k (l + hb k ) - a fe _i(l + hb k -\), a k = a k (b k+1 - b k + ha k+1 - ha k ) (3.51) 



2. In the coordinates (a, 6) the linear Poisson bracket ( |3.43p of the RTL formally coincides 
(incidentally) with the linear Poisson bracket ( |3.2|) of the usual Toda lattice: 



{b k , a k }i — —a k , {a k ,b k+1 }i 



a k 



(3.52) 



Note also that in the coordinates (a, d) the quadratic Poisson structure ( |3.44| ) formally 
coincides with the quadratic Poisson bracket (|3.4j) for the usual Toda lattice. 



3. The third ("cubic") invariant Poisson structure is known for the RTL hierarchy, how- 
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ever its r-matrix interpretation is unknown: 



{a>k, a kj3 — 


—a k [a k -t a k -\- nacb) 


{flfc,<ifc+i}3 = 


-a k (d 2 k+1 + a k + hd k+1 a) 


{dk,d k+ i} 3 = 


— a k(d k + d k+ i + hd k d k+1 ) 


{a>ki «fc+l}3 = 


— o, k a k+ \{2d k+ i + /id/. + ha k+ \ 


{rffc, 0-^+1)3 = 


-a k a k+1 (l + /id*) 


{ a k^d k+2 } 3 — 


-o fc o fc+ i(l + /id fc+2 ) 


\0,k, 0, k+2 } 3 = 


—ha k a k+ \a k+ 2 



(3.53) 



4 Volterra and relativistic Volterra lattices 
4.1 VL 

We consider here the equations of motion of VL (|1.2| ) under periodic boundary conditions, 
so that the corresponding phase space is: 

V = R 2N (u 1 ,v 1 ,...,u N ,v N ) (4.1) 

There exist two compatible local Poisson brackets on V invariant under the flow VL: 

{u k ,v k } 2 = -UkVk, {v k ,u k+1 } 2 = -v k u k+1 (4.2) 

and 

{u k} V k } 3 = -U k V k {u k + V k ), {v kl U k+1 } 3 = -V k U k+ i(v k + Uk+l) 

{u k ,u k+1 } 3 = -u k v k u k+1 , {vk,Vk+ih = -v k u k+1 v k+1 

The corresponding Hamilton functions for the flow VL are equal to 

N N 

K 1 (u,v) = ^u k + ^v k (4.4) 

k=l k=l 

and 

N N 

Bo(u,v) = ^log(M fc ) or R (u,v) = J2 l °s(vk) (4.5) 

k=l k=l 
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(4.3) 



(the difference of the latter two functions is a Casimir of {•, -}a). 

The Lax representation of VL we use here lives in g = g ® g and was introduced in 
|[K1|| (where the system ( |1.2j ) was called "modified Toda lattice"), see also |55| . Consider the 
following two matrices: 

TV N 

U (u, v, A) = A Ek+i,k + u kEk,k = £ + u (4.6) 

k=l k=l 

N 

V{u, v, A) = I + A" 1 £ v k E k>k+1 = 1 + vE- 1 (4.7) 

k=l 

These formulas define the "Lax matrix" (U, V) : V i— > g <S> g- The flow ( p.. 2D is equivalent to 
either of the following Lax equations in g ® g: 

f U = UB 2 -B 1 U f U = A 1 U-UA 2 

\ V = VB X - B 2 V [V = A 2 V -VAx 1 ' ] 

with the matrices 

Bl = n+ (UV), B 2 = n + (VU), A x — tt^(UV), A 2 = tc.(VU) (4.9) 

so that 

N N 

B l = A ^2 E k+ i ik + ^2(u k + v k _i)E kk (4.10) 

k=l k=l 

N N 

B 2 = \Y J E k+i,k + Y.^ u k + v k)E k k (4.11) 

k=l k=l 

N 

A x = A -1 u k v k E k)k+x (4.12) 

k=l 
N 

A 2 = \~ l ^u k+1 v k E k , k+1 (4.13) 

k=l 

As a consequence, the matrices 

T 1 (u,v,X) = U(u,v,\)V{u,v,\), T 2 (u,v,X) — V(u,v, X)U(u,v, X) (4.14) 
satisfy the usual Lax equations in g: 

f=[T,±n±(T)] (4.15) 



The Lax equations ( |4.8| ), ( |4.15| ) may be given an r-matrix interpretation. 
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Theorem 4.1 ||S5|| The Lax matrix map (U, V) : V ^ g ® g is Poisson, if V is equipped 
with the bracket {•, -} 2 and g (g> g is equipped with the bracket PB 2 (A l5 A 2 , S) corresponding 
to the operators A 1; A 2 , S defined in ( ft.49| ). The maps T 12 : g ® g ^ g, 



T 1 (U,V) = UV, T 2 (U,V) = VU 
are Poisson, if the target space g is equipped with the bracket PB 2 (A 1; A 2 ,S) 

It is easy to see that the matrices Ti,T 2 formally coincide with the Lax matrix (|3.8| ) of the 
To da lattice, with the coordinates (a, b) given by the corresponding formula for the Miura 
maps in ( |1.8|) . The Poisson property of the monodromy map is therefore reflected in the 
(first half) of the following statements about the Miura maps. 

Theorem 4.2 The Miura maps -Mi j2 : V i— > T are Poisson, ifV carries the bracket {•, -} 2 
( |4.2|) and T carries the bracket {•, -} 2 ( p. 41) , and also ifV carries the bracket {■, ( [4.3|) and 
T carries the bracket {-, -} 3 ( |3.6|) . 



4.2 VL -» dVL 

To find an integrable time discretization for the flow VL, we apply the recipe of Sect 
with F(T) — I + hT, i.e. we consider the map described by the discrete time "Lax triads" 

(U = B?UB 2 (U = Al UA^ 

\ V = B 2 1 VB l ° r \V = A 2 VA^ l4 ' ibj 

with 

B 1 = IL+(I + hUV), B 2 = U + (I + hVU), A 1 =U-(I + hUV), A 2 = Ii^{I + hVU) 



Theorem 4.3 Consider the change of variables V(u, v) i— > V(u, v) defined by the formulas 

u k = ujb(l + hvk-i), v k = v k (l + hu k ) (4.17) 

The discrete time Lax equations ( [4.1 6|) are equivalent to the map (u, v) \— > (w, which in 
coordinates (u, v) is described by the following equations: 

u fc (l + /ivfc-i) = u fc (l + kv k ), vjb(l + Mi fc ) = v fc (l + /iu fc+ i) (4.18) 
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Proof. The formulas (|4 . 1 7 ) allow to find the factors -B 12 , A 12 in a closed form. Indeed, 



with the help of ( |4. 1 7| ) we can represent ( |4.6| ), ( |4.71 ) as 

TV TV 

U = J2u k {l + hvk^Ekk + XY.Ek+i^ (4.19) 

k=l k=l 

TV 

V = / + A- 1 ^v fe (l + / i u fe )£; ifc)fe+1 (4.20) 

k=l 

From these formulas we derive: 

TV TV 

I + hUV = h\J2 E k+i,k + J2 ((l + ^u fe )(l + /iv fc _i) + h 2 u k - 1 v k - 1 )E kk 

fc=l k=l 

TV 

+h\~ 1 u fc Vfe(l + hu k )(l + /iv fc _i)E fejfc+ i 
fc=i 

Obviously, this matrix may be factorized as / + hUV = BiA±, where 

TV TV 

B l = Y.( l + h "k)(l + hv k _ l )E kk + h\Y J E k+1 , k (4.21) 

k=l k=l 

TV 

A l = I + hX' 1 u k \r k E ktk+1 (4.22) 

k=l 

Similarly, we find: 

TV TV 



I + hVU = /iA^^fc+i,ifc + Z)((l + ^u fe )(l + /iv fc ) + /i 2 u fe v fe _i) J E;, 

k=l k=l 

TV 

+h\~ 1 Y u fc+ iv fe (l + hu k )(l + hw k )E k}k+l 



kk 



k=l 

which may be factorized as I + hVU = B2A2, where 

TV TV 

B 2 = J2( 1 + hu k)(l + hv k )E kk + h\J2E k+ltk (4.23) 

k=l k=l 

TV 

A 2 = I + hX-^Uk+iVkE,^ (4.24) 

k=l 

Now the discrete time Lax equations B-JJ = UB2, B 2 V = VB\ (or their A-analogs) 
immediately imply the equations of motion ( f4.18|) . I 
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Comparing two pairs of formulas (|4.21|) , (|4.22j ) and Q4.23Q , (|4.24j ) for dVL with the 
formulas fl3.15| ), ( |3.16| ), we see immediately that the Miura maps Aii^ '■ V(w, v ) i— > T(a, 6) are 
conjugated by the changes of variables ( |4.17|) , ( |3.13| ) with the maps M 12 : V(u, v) i— > T(a, b) 
given by: 



Mi : 



i + hh k = (i + /iUfe)(i + /iv fe _ 1 ; 

a k = u k v k 



Mo 



l + hb k = (l + hu k )(l + hv k ) 



(4-25) 

The map ( f4.18|) , denoted hereafter by dVL, serves as a difference approximation to the 
flow VL (|1.2j ). By construction, in the variables (u, v) this map shares the integrals of motion 
and the invariant Poisson structures with the flow VL. In the coordinates (u, v) the integrals 
of motion and the Poisson brackets become deformed, the latter ones become nonlocal, since 
the inverse of the map ( |4.17| ) is nonlocal. Nevertheless, we have the following local invariant 
Poisson bracket for the map dVL: 

Theorem 4.4 | |K1[ ] The pull-back of the bracket 

{-,-} 2 + M-,-}3 (4.26) 
on V(u, v) under the change of variables ( |4.17|) is the following bracket on V(u, v): 
{u fc ,v fc } = -u fc v A .(l + /m fc )(l + /iv fe ), {v fc ,u fc+ i} = -v fc u fe+1 (l + /iv fc )(l + /m fc+ i) (4.27) 

The map ( |4.18| ) is Poisson with respect to the bracket ( f4.27|) . The maps M li2 (|4.25|) are 
Poisson, ifV(u,v) carries the bracket ( f4.27| ) and T(a, b) carries the bracket ( |3.21| ). 

Proof - by a direct verification. ■ 

The map (|4.18| ) was found in [THO|| , |[HTI |, where, however, its place in the continuous 
time Volterra hierarchy and its Poisson structure were not elaborated. The Miura maps 
( |i.25| ) were found in JKI] without connecting them to the discretization problem, and also 
in ETj in the discretization context, but without Poisson properties. 

The previous theorem allows also to calculate the equations of motion of continuous time 
hierarchy, to which the map dVL belongs. 

Theorem 4.5 |[K1|| The pull-back of the flow VL under the map (|4.17|) is described by the 
following equations of motion: 



u k = u fe (l + hu k )(vk ~ Vfc-i), v fc = v fc (l + hv k )(u k+1 - vjfc) 



(4.28) 
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Proof. Since the flow VL has a Hamilton function h~~ l logHo(w, v) in the Poisson bracket 
h{-, and this function is a Casimir of the bracket {•, -}2, we conclude that this function 
generates the flow VL also in the bracket ( |4.26|) . This means that the pull-back of the flow 
VL is a Hamiltonian flow in the bracket (|4.27|) with the Hamilton function 



N 



h 1 ^ log (u fc (l + /tv fe _ 



k=l 



Calculating the corresponding equations of motion, we arrive at the equations (|4.28| ) 



4.3 dVL explicit dVL 

We use the same trick as in the case of dTL, in order to extract from dVL the explicit 
discretization for VL. Namely, we attach the variables = u&(t), = Vfc(i) to the lattice 
site (k,t) G Z x hZ, and then perform the change of independent variables (k,t) i— > (k,r) = 
(k, t + kh). We denote 

itfe(r) = ujfc(t) — u k (r — kh), v k (r) = v k (t) = v fc (r - kh) (4.29) 

Denoting the h— shift in r still by the tilde, we immediately derive from ( f4.18| ) the following 
difference equations for the variables u k , 

G fc (l + hv k _ x ) = u fc (l + hv k ), v k (l + hu k ) = v k (l + hu h+x ) (4.30) 

This is clearly an explicit discretization, since, knowing (it, v), it allows to calculate explicitly 
first u, and then v. 

We find now Lax representations for the map (|4.30|) . As in the case of the dTL, there 
exist two versions thereof. 



Theorem 4.6 The map ( 4.30|) allows the following Lax representation: 



U = C 1 UC 2 , VC l = C^ l V (4.31) 



with the matrices 



N N 

U(u,v,X) = Y, u kEkk + Xj2 E k+i,k (4.32) 

k=l k=l 
N 

V(u,v,X) = J + A-^vjfe^fc+i ( 4 -33) 

k=l 
N 

C{u,v,X) = I - hX' 1 u k v k E ktk+1 (4.34) 

k=i 

N 

C 2 (u, v, X) = I - hX' 1 u k+1 v k E ktk+1 (4.35) 

k=i 



27 



The Lax representation ( }4.31|) implies also 

L = C l L C UVC' 1 = C X UV 

with the Lax matrix 



(4.36) 
(4.37) 



L = UVC 1 

Proof. We start with the following Lax representation of the dVL in g <g> g: 

U = A^A- 1 , V = A 2 VA- X 

with the matrices ( |4.19| ), ( |4.20| ) , ( |4.22| ), fl4.24| ). These Lax equations may be considered as 
the compatibility condition of the following linear problems: 

U(f> = fiip { 4> — A 2 (, 



ip = Attp 



We shall call the first pair the spectral linear problem, and the second pair the evolutionary 
linear problem. In coordinates the above equations may be presented as 

h = 4>k + /iA _1 u fe+ iv fe A . +1 



U fc (l + /lV fc _i)</>& + A0 



k-1 



k + A ^(1 + hu k )i) k+1 = fj4 k { ip k = < l p k J r h\ 1 u fc v A .<^ fc+1 

Upon the change of variables (k, t) i— > (k, r) this implies: 

f u k (l + h v k _ x )<j> k + A c? k _ x = ^ k ( ^ = ^ + hX -iu k+1 v k $ k+1 



^ Vfe + A l v k (l + hu k )il> k+l = [i<f> k 



^k = ^k + ^A 1 u k v k cf) k+1 



After simple manipulations (substitute into the first pair of equalities the values of 4> k _i and 
i/) k obtained from the second pair of equalities) this can be brought into the form 
' u k cj) k + A0 fc _! = fiip k { cf) k = 4> k - h\~ 1; u k+1 v k <i> k+1 

V'fc = i>k - h\- l u k v k 4) k+l 



k ip k + A 1 v k tp k+1 = fi(j) k 
or in the matrix form 



U<fi = flip 



where we introduced the matrices U, V, C, C 2 from ( f4.32|) -( |4.35|) . The compatibility 
condition of these linear problems coincides with ( [4. 31 ). ■ 
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Theorem 4.7 The map (|4.30|) allows the following Lax representation: 

f- l UA 2 = (A^Ag l T) ■ T~ X UA 2 ■ (^AxAa) 

A 2 l V = (A 2 'A^f) • A^V • GF^AgAO 

wt/i i/ie matrices U , V , T from ( [4. 32|) , Q4.33|) , ( |3.35|) , and 

Ai = diag(l + /iu fc ), A 2 = diag(l + /iv fe ) 

A 3 = diag(l + fru fc _i) 
Tne Lax representation ( |4.38|) implies also 

L + = (A^A- 1 ?) ■ L + ■ (J" _1 A 3 Ai) 

mt/j t/ie Lax matrix 

L + = F~ l UV 

Proof. This time we start with the following Lax representation of the map dVL: 

U = Bi 1 UB 2 , V = B 2 1 VB 1 



(4.38) 



(4.39) 
(4.40) 

(4.41) 
(4.42) 



with the matrices ( |4.19| ), ( |4.20| ), ( |4.21|) , (f4.23|) . These Lax equations are the compatibility 
conditions of the following linear problems: 



U(p = nip 



B 9 



(it turns out to be more convenient to take the tilded versions of the equations in the first 
pair). In components the above equations take the form: 

u fe (l + /iv fc _i)0 fe + A0 fc _i = fi^ k ( (1 + hu k )(l + hv k )4> k + /iA0 fc _i = <p k 

i) k + A -1 v fc (l + h\i k )i} k+1 = n4> k { (1 + frv fe _i)(l + h\i k )4) k + hX^i = ip k 

After the change of variables (k, t) \— > (k, r) this takes the form: 

u k (l + hv k _ 1 )(f> k + \(/> k _ 1 = /jif> k I (l + hu k )(l + hv k )(f> k + h\(l) k _ 1 = (j> k 



i\) k + A l v k (l + hu k )4> 



k+l 



[1 + h Ufe-i)(l + hu k )tp k + hXif)^ = ij> h 



29 



Now we use the equations of motion ( f4.30|) to bring the first pair (the spectral problem) into 
the form ^ _ _ 

u k (l + hv k )4> k + A0 fc _ 1 = /iV'fc 

(4.43) 

i> k + + hu k+1 )^ k+l = [i(f> k 

The second pair (the evolutionary problem) may be presented as 



1 + hu k )(l + hv k )<f> k = tf> h - hXcf)^ 
[1 +hv k - 1 )(l + hu k )^ k = if) k - hXif)^ 



<t> = A 2 " 1 Ar 1 ^ 



1p = A^A^Fl/) 



(4.44) 



Now we transform the equations describing the spectral problems. From the first equations 
of the two pairs in ( 4.43| ), ( 4.44Q it is easy to derive: 



0fe - (1 + hv k)4> k = h ^k 
and using this in the first equation in ( 4.43|) , we bring the latter into the form 

u k (l + hv k )4> k + h\{\ + hv k -i)4> k -i = /i(V? fc - h\^) k _{) 
From the second equation in ( [4.44 ) we find the value 

1 



(4.45) 



(1 + hu k+1 )ip k+1 



-Ofc+i - h\i/> k ) 



1 + hv k 

which, being substituted into the second spectral equation in ( 4.43|) , implies: 

%h + X ~ lv k^ k +i) = V<f> k 



1 + hv k 

Now we can put the spectral problems in ( 4.45|) , (|4.46| ) into the matrix form: 

F^UA-rf) = 1$ 

A^Vj; = ^ 

The compatibilty conditions of the linear problems (|4.44|) , (|4.47|) coincide with (|4.41 ). 



(4.46) 



(4.47) 
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4.4 Explicit dVL -> RVL 

Now we have obtained two Lax matrices for the explicit dVL: 

L-(u, v, A) = UVC' 1 and L+(u, v, A) = F^UV 
The relation between them is the following: 

Q l U{X)Q 2 1 = V T {\'), Q 2 V{X)Qi 1 = £/ T (A') ftiC(A) ftf 1 = ^(A') 

so that 

fiiZ_(A) fif 1 = Lj(A') 
where A' = a -1 A -1 , a = {uiVi...u N v N ) 1 / N and 

Qi = diag(l, a~ l u-iVi, a~ 2 u 1 v 1 u 2 v 2 , . . . , a~ N+1 UiV i...%_i%_i) 
fi 2 = diag(l, a _1 viti2, aT 2, Uiit 2 v 2 M3, • • • , a _Ar+1 Viit 2 ---f Af-i^Ar) 

We postulate that the hypothetic relativistic Volterra lattice (RVL) lives in the hierarchy as- 
sociated to these Lax matrices, and it remains to identify this hierarchy. (The above relation 
between L± assures that the both Lax matrices generate one and the same hierarchy). We 
start with the Lax matrix L . 

Theorem 4.8 The triples (U, C^ 1 , V) form a Poisson Lax matrix map V(u, v) i— > g®g®g, 
if V(u, v) carries the Poisson bracket 

{u k , v k } 2 = -u k v k , {v k , u k+ i} 2 = -v k u k+l (4.48) 

and g ® g ® g carries the Poisson bracket PB 2 (Ai, A 2 , S) defined by the operators 

/ Ai -S -S \ / A 2 -S* -S* \ / S S S \ 

Ai = S* Ai S* , A 2 = S A 2 -S* , S = S -S* -S* 
V s* -S A x J V s s A 2 / V s s -S* / 

(4.49) 

The Lax matrix maps L = UVC^ 1 , VC~ 1 U ', C~ 1 UV : V(u, v) i— > g are Poisson, if the 
target g is equipped with PB 2 (Ai, A 2 , S). The hierarchy of continuous time flows is described 
by the Lax triads 

U = ±U-7i ± (f(VC- 1 U))Tir ± (f(UVC- 1 ))-U (4.50) 
C = ±C-7r ± (f(C- l UV))Tn ± (f(UVC- l ))-C (4.51) 
V = ±V-n ± (f(C- l UV))TK±(f(VC~ l U))-V (4.52) 
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As a consequence, the evolution of the matrix L = UVC 1 is described by the standard 
Lax equation 

L = [X_,±7r±(/(X_))" 

In particular, the "first" flow corresponding to f(L) = L, coincides with the RVL ( |1.13| ), 
and the auxiliary matrices in the Lax triads are given by 

N N 

n+iUVC- 1 ) = J2( u k + v k-i + hu k _ 1 v k _ 1 )E kk + \J2E k+1 , k (4.53) 

fc=l k=l 
N N 

n+iC^UV) = Y,( u k + v k-i + hu k v k )E kk + Xj2 E k+i,k (4.54) 

k=l k=l 
N N 

7T + {VC- l U) = Y,( u k + v k + hu k v k )E kk + \Y,E k+1>k (4.55) 

k=l k=l 

The map ( |4.30| ) belongs to this hierarchy (in particular, is Poisson with respect to the bracket 
( |4.48| ) ) and corresponds to f(L) = h^ 1 log(I + hL). 

Proof. The first statement is proved with the help of the tensor notations for the quadratic 
r-matrix brackets, along the same lines as the proof of analogous statements in ||S2| , |[55|| . 
All other statements, except the last one, are consequences of the first one and Theorem El) 



To identify the place of the map (|4.30|) in this hierarchy, we have, referring to Theorem ^L6 
to prove the following relations: 

cr 1 = n_(i + huvc- 1 ), c 2 x = n_(/ + hvc~ x u), c' 1 = + hc^uv) 

Clearly, it is enough to prove only the first of them, but it follows immediately from 

/ + hUVC- 1 = (C+ hUV)C- x 

and 

N N 

C + hUV = J2(± + hu k + hv k ^)E kk + h\J2 E k+1>k e G + , C- 1 e G 

k=l k=l 

The proof is complete. ■ 

Turning to the case of the Lax matrix L + , we have the following results. 

Theorem 4.9 The hierarchy associated with the Lax matrix 

L, = T' X UV 
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is described by the following Lax triads: 

If = AxU - UA 2 , V = A 2 V - VA (4.56) 
which have to be supplemented by the identity 

F = = A\T - TA Q (4.57) 
so that for the matrix L + = J r ~ 1 UV there holds the Lax equation 

L + =[Aq,L+] (4.58) 

Here 

A = n_(f(L + ))-o-(f(L + )) (4.59) 

and the matrices A\$ are uniquely defined by the condition of compatibility of the equations 
of the Lax triads (|4.56|) , (|4.57|) . In particular, the "first" flow of the hierarchy corresponding 
to f(L) = L, coincides with RVL (|1.13| ), and for this flow 

N N 

A = A -1 UkV k Ek,k+i ~hJ2 u k v k E kk (4.60) 

k=l k=l 
N N 

•Ai = A -1 u k v k E ktk+1 - h u k _ x v k _ x E kk (4.61) 
fc=i fc=i 

N N 

M = A -1 u k+1 v k E Kk+ i -hJ2 u k v k E kk (4.62) 
fc=i fc=i 

Proof goes by a direct check. We define Aq by ([4.59| ), then put A\ = TA§T~ X G g © g . 
The matrix A 2 G g © g may be defined in two ways to make either of two Lax triads 
in Q4.56| ) consistent. These two definitions have to be consistent, because the choice ( 4.59Q 
assures the validity of the Lax equation ( ^4.58| ) for the matrix L + , due to Theorem |3.8[ I 



Comparing the Lax matrices L± for the RTL hierarchy (see ( |3.40| , ( |3.50D ) and for the RVL 
hierarchy (see ( f4.36| ), ( 4.41|) ), we see that the Miura map between two hierarchies formally 
coincides with the Miura map between the TL hierarchy and the VL hierarchy, because it is 
equivalent to identifying L(a, d, A) with U (u, v, A) V(w, v, A), where L, U, V are the matrices 
from the nonrelativistic theories. This Miura map is 



Mi(u, v) = (a, d) : 



d k = u k + v k _i 
a k = u k v k 



This Miura map is Poisson with respect to the invariant quadratic Poisson brackets, since 
these also formally coincide with the invariant quadratic Poisson brackets of the nonrela- 
tivistic hierarchies. Recall that also the cubic Poisson bracket (|3.53|) is known for the RTL 
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hierarchy. It turns out that its pull-back by the Miura map is given by the following nice 
formulas: 

{u k , v k } 3 = -u k v k (u k + v k + hu k v k ), {v k , u k+1 } 3 = -v k u k+ i{v k + U k+ i) 

{u k , u k+1 } 3 = -u k v k u k+1 (l + hu k ), {v k , v k+1 } 3 = -v k u k+1 v k+1 (l + hv k+1 ) 
{v k ,u k+2 } 3 = -hv k u k+1 v k+1 u k+2 (4.63) 



Finally, let us comment on the twin RVL flow (|1.14j) . Obviously, this system, as well as the 



corresponding Miura map A4 2 , m ay be obtained upon the renaming u k — > v k , v k — > u k +i. 
The Lax matrices for the corresponding hierarchy arise from the ones previously concerned 
upon change UV to VU. The explicit discretization of the VL living in this hierarchy reads: 

u k (l + hv k -x) = u k (l + hv k ), v k (l + hu k ) = v k (l + hu k+ i) (4.64) 

(here, knowing u, v, one calculates first v, and then u). 

5 Bogoyavlensky lattices 
5.1 BL2 

A class of systems called Bogoyavlensky lattices (although discovered several times indepen- 
dently before his papers |B|, see the references in and also ||K1|| ) serve as generalizations 



of the Volterra lattice. We consider here only one representative of this class, which is 
however typical. 

The integrable lattice system considered in the present section will be called BL2. Its 
phase space (in the periodic case) is 

B = R 3N (u 1 ,v 1 ,w 1 , ...,u N , v N , w N ) (5.1) 

The equations of motion read: 

Uk = u k (v k + w k - V k - X - W k -l) 

v k = v k (u k+ i + w k -u k - w k -i) (5.2) 
w k = w k (u k+1 + v k+1 - u k - v k ) 

This system is Hamiltonian with the local Poisson brackets on B 

{u k ,v k } 2 = -u k v k , {v k ,u k+1 } 2 = -v k u k+l 
{u k ,w k } 2 = -u k w k , {w k ,u k+l } 2 = -w k u k+ i (5.3) 
{vk,w k } 2 = -v k w k , {w k ,v k+1 } 2 = -w k v k+1 
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and the Hamilton function 



(5.4) 



N N N 

Kl(u, V, W) = U k + V k + H W k 
k=l k=l k=\ 

Notice that the system BL2 is more familiar in the form 
which arises upon the re-naming of variables 

Uk = Cl3k-2, Vk = 03fe-l, W k = d 3k 

We use here the Lax representation of BL2 living in g = g <g> g <g> g and introduced in 
~5l. Consider the following three matrices: 



U (u, v, w, A) 
V(u, v, w, A) 
W(u,v,w, A) 



N N 

A Y E k+l,k + Y U k E k,k = £ + u 
k=l k=l 

N 

I + X' 1 ^ v k E k>k+1 = 1 + vS- 1 

k=l 
N 

I + A" 1 Y w kE k ,k+i = I + w£~ l 
k=i 



(5.5) 
(5.6) 
(5.7) 



These formulas define the "Lax matrix" (U,V,W) : B 1— > g <8> g <E> g- The flow (|5.2|) is 
equivalent to either of the following Lax equations in g <g> g <g> g: 



U = UB 3 - B X U 
V = VB 1 - B 2 V 
W = WB 2 - B 3 W 



or 



with the matrices 

B x = n+(UWV) , B 2 = n+(VUW) , 
A x = vr_ (UWV) , A 2 = vr_ (VUW) , 



U = A X XJ - UA 3 
V = A 2 V - VA X 
W = A 3 V - VA 2 



B 3 = ir + (WVU) 
A 3 = tt_(WVU) 



(5.8) 



so that 



JV N 

B\ = A Y E k+ i :k + J2( u k + v k-i + w k -i)E kk 

k=l k=l 
N N 
B 2 = \J2 E k+l,k + J2( U k + V k + 'W k -l)E kk 
k=l k=l 



(5.9) 
(5.10) 

(5.11) 
(5.12) 
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N N 



B 3 = A E k+i,k + J2(u k + v k + w k )E kk (5.13) 

k=l k=l 

N N 

Ai = X~ 1 ^2(w k ^xv k + u k v k + u k w k )E kjk+1 + \~ 2 ^2ukW k v k+1 E kik+2 
k=i fe=i 

N N 

A 2 = A -1 J2(u k w k + v k w k + v k u k+ i)E Kk+l + A~ 2 ^ v k u k+1 w k+1 E k:k+2 

k=l k=l 
N N 

A?, = A -1 J2(v k u k+1 + w k u k+1 + w k v k+ i)E Kk+1 + A~ 2 ^ w k v k+1 u k+2 E kik+2 

k=l k=l 



The Lax equations (|5.8| ) allow the following r-matrix interpretation. 

Theorem 5.1 ||S5|| The Lax matrix map (U, V, W) : B \— > g®g®g is Poisson, ifB is equipped 
with the bracket {•, -} 2 and g®g®g is equipped with the bracket PB2(Ai, A2, S) corresponding 
to the operators A 1; A 2 , S defined in ( |4.49| ). The monodromy maps T 1)2j3 : g ® g ® g 1— > 

Ti (£7, V, W) = £7 W, T 2 (£7, V, W) = VE/W, T 3 ([/, V, W) = W VU 

are Poisson, if the target space g is equipped with the bracket PB 2 (Ai, A 2 , S). 

5.2 BL2 -> dBL2 

To find an integrable time discretization for the flow BL2, we again apply the recipe of Sect. 
|2.3| with F(T) — I + hT, i.e. we consider the map described by the discrete time "Lax 
triads" 

' U = B^UB 3 { U = AxUA^ 1 

(5.14) 



u = 


B^ l UB 3 






All/As 1 


V = 


B- X VB X 


or < 




A 2 VA- 1 


w = 


-- B^WB 2 




k w = 


-- A 3 WA^ 



with 

Bx = U+(I + hUWV) , B 2 = IT+ (7 + WE7W) , B 3 = n+ ( J + h W£7) 
Ax = II_(I + /iE7W), A 2 = U_(I + hVUW), A 3 = n_(7 + hWVU) 



Theorem 5.2 Consider the change of variables B(u, v, w) 1— > B(u,v,w) defined by the for- 
mulas 



U k = U fe (l + /lWjfe_i)(l + /iv fc _i) 

f fc = v fc (l + /m fc )(l + hw k _i) (5.15) 

Wfc = Wfc(l +/iV fe )(l + /lUjfc) 
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The discrete time Lax equations Q5.14Q are equivalent to the map (w, v, w) h- > (u, v, w) which 
in coordinates (u, v, w) is described by the following equations: 

u k (l + /ivjfe_i)(l + hwk-i) = u fc (l + hv k )(l + hw k ) 
v fc (l + ft,Ufc)(l + hw k -i) = v fc (l + ft,u fc+ i)(l + /iwfe) (5.16) 
w fc (l + fru fe )(l + /iVfc) = w fc (l + hu k+1 )(l + /iVfc+i) 

Proof is parallel to the proof of Theorem [4 . 3| . The formulas (|5.15 ) allow to find the factors 
.61,2,3, ^-1,2,3 in a closed form. Indeed, with the help of ( |5.15| ) we can represent (|5.5|)-([577D 

as 

TV TV 

U = ^u fc (l + /iw fc _ 1 )(l + /iv fc _ 1 )£ fcfc + \J2 E k+i,k (5.17) 

k=l k=l 

TV 

V = J + A- 1 $> fc (l + /m fc )(l + hw k ^)E ktk+ i (5.18) 

k=l 

TV 

W = I + \- ly £w k (l + hv k )(l + hu k )E kthH (5.19) 

k=l 

From these formulas we derive the following expressions for the factors Bi = n + (7 + hTj) 
and Ai = II_(J + /iT<), £ = 1,2, 3: 

AT TV 

= ^(l + /iu fc )(l + /iw fc _ 1 )(l + /iv fc _ 1 ) J E fefc + /iA^E fc+ljfc (5.20) 
fc=i fe=i 

TV TV 

£ 2 = ^(l + / i v fc )(l + ^u fc )(l + /iw fc _ 1 ) J E fefc + / i A^E fc+lifc (5.21) 

k=l k=l 

TV TV 

B 3 = ^2(l + hw k )(l + hY k )(l + hu k )E kk + h\J2E k+ i, k (5.22) 

fc=l k=l 
TV 

Ai = / + hX' 1 ^(w fe _iv fc + u fe v fc + u fc w fe + /iw^iUfcVfc + hu k v k w k )E ktk+1 

k=l 

TV 

+/iA -2 ^ u fe w fe v fe+1 (l + /iv fc )(l + /iw fc )(l + hu k+l )E Kk+2 (5.23) 
fe=i 

TV 

A 2 = I + hX' 1 ^2(u k w k + v k w k + v k u k+1 + hu k v k w k + hv k w k u k+1 )E k , k+1 

k=l 

TV 

+hX~ 2 v fc u fe+1 w fe+1 (l + hw k )(l + hu k+1 )(l + hv k +i)Ek,k+2 (5-24) 
fe=i 
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N 



A 3 = I + hX 1 J2(v k u k+1 + w fc u fc+ i + w fc v fe+ i + hv k w k u k+1 + hw k u k+1 v k+1 ) E kik+1 

k=l 

N 

+h\~ 2 w fc v fc+ iu fc+2 (l + hu k+1 )(l + kv k +i)0- + hw k+ i)E k)k+2 (5.25) 



k=l 



Now the equations of motion ( 5.16Q follow directly from the B version of the discrete time 
Lax triads (|5.14p. ■ 



The map (|5.16|) , denoted hereafter by dBL2, serves as a difference approximation to the 
flow BL2 ( |5.2j) . According to the general properties of discretizations of Sect. [2.3| , this map 
in the coordinates (u, v,w) is Poisson with respect to the bracket (|5.3|). Unfortunately, this 
bracket becomes highly nonlocal in the coordinates (u, v, w). After the re-naming 

Ufc = a 3 fc_ 2 , Vfe = a 3 fc_i, Wfe = a 3 fc 

the map dBL2 turns into 

a*(l + ha*-i)(l + /ia fc _ 2 ) = a*(l + ha k+1 )(l + ha k+2 ) 

and in this form it appeared for the first time in JTHD|, and then in [[PN| , [53|. 



Wfc(r) 


= Ufc(t) 


= u fc (r - 


kh) 


v k (r) 


= v fc (t) 


= v fc (r - 


kh) 


w k (r) 


= w fc (t) 


= w fe (r - 


- kh) 



5.3 dBL2 -> explicit dBL2 

Following the scheme used already for dTL and dVL, we extract an explicit discretization 
for the flow BL2 from the map dBL2. To this end we denote 



(5.26) 



These variables satisfy the following difference equations with respect to the discrete time r: 

u k (l + hv k -i)(l + hw k _i) = u k (l + hv k )(l + hw k ) 
v k (l+hu k )(l + hw k -i) = v k (l + hu k+1 )(l + hw k ) (5.27) 
w k (l + hu k )(l + hv k ) = w k (l + hu k+1 )(l + hv k+1 ) (5.28) 

This is indeed an explicit discretization, because, knowing (n, v, w), one calculates explicitly 
u, v, w (in this order). We now turn to finding a Lax representation for this map. Unlike 
the previous cases, we could find reasonable formulas only starting from the .B-version of 
the Lax representation for dBL2. 
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Theorem 5.3 The map (|5.27|) allows the following Lax representation: 

{ J^- 1 UA 2 A 3 = (A^A^A^J 7 ) ■ f- 1 UA 2 A 3 ■ (^A^As) 

A 2 X V = (A^A^Ag 1 T) ■ A^ X V ■ (^- 1 A 5 A 4 A 1 ) 

A^W = (A^A^A^J 7 ) • A3 1 W • (JT- 1 A 5 AiA 2 ) 



wi/i £/ie matrices 



N N 

U = X ^ Ek+l,k + X] U kEk,k 
k=l k=l 



N 



V = I + X 1 Y, v kE k ,k+i 

k=l 
N 

W = I + \- 1 Y.w k Ek, k +i 



k=l 



and 



Ai = diag(l + huk), A 2 = diag(l + hv k ), A 3 = diag(l + hw h 

A 4 = diag(l + hv k -i), A 5 = diag(l + hw k -i) 
The Lax representation ( 5.29|) implies also 

L + = (A^A^As 1 ^) • L + ■ (^- 1 A 5 A 4 A 1 ) 

with the Lax matrix 

L + = (J r - 1 UA 2 A 3 )(A 3 1 W)(A 2 1 V) = J r - 1 UA 2 WA 2 1 V 
Proof. We start with the following Lax representation of the map dBL2: 
U = B^UB 3 , V = B 2 1 VB 1 , W = B^WB 2 



(5.29) 



(5.30) 
(5.31) 
(5.32) 

(5.33) 
(5.34) 

(5.35) 
(5.36) 



with the matrices ( |5.17|) -( [5.19|) , ( |5.20| )- (|5.22| ). These Lax equations are the compatibility 
conditions of the following linear problems: 



Ucj) = flip 



f B ? 



< b ± 4j = ip 

. B 2 x = X 
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(it turns out to be more convenient to take the tilded version of the spectral problem, i.e. of 
the equations in the first triple). In components the equations of the spectral linear problem 
take the form: 

Ufc(l + /lWfc_i)(l + /lV fc _i)0 fc + A0fe_i = pLIpk 

i) k + A _1 v fe (l + hu k )(\ + /iwn)^ + i = nxk 

, Xk + A _1 w fc (l + hv k )(l + hu k )xk+i = lM>k 
while the equations of the evolutionary linear problem take the form 

(1 + hu k )(l + hv k )(l + hw k )cj) k + /iA0 fc _i = <p k 

[1 + /iv fc _i)(l + /iw fc _i)(l + hM k )i) k + /iA^ fc _i = ip k 

k (1 + /iw fc _i)(l + /iu fc )(l + /iv fc )xfc + hAxfc-i = Xfc 
After the change of variables (fc, t) h- > (/c, r) the spectral equations take the form: 

tt fc (l + /iw fc _i)(l + hw k ^)4) k + A0 fe _ 1 = 

^ fc + A _1 « fc (l + + /i™ fe _i)V> fc+1 = 

k X k + A _1 w? fe (l + /i« fc )(l + hu k )x k+1 = fi<t> k 
and the evolutionary equations take the form 

(1 + hu k )(l + hv k )(l + hw k )(p k + /iA0 fc _! = fe 



1 + /i« fc _x)(l + /iw M )(l + hu k )if) k + hXif)^ = *l>k 



[1 + hw k -i)(l + huk)(l + hv k )Xk + h\Xk-i = Xk 



(5.37) 



or, in the matrix notations, 



= A^A^A^Ff 
j> = Af 1 A^Af Fij> 

X = A 2 1 A- l 1 A^F X 



(5.38) 
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Now we use the equations of motion (|5.27|) to bring the spectral equations into the form 

u k (l + hv k )(l + hw k )4> k + A0 fc _ 1 = (iij) h 



tp k + X 1 v k (l + hw k )(l + hu k+1 )ip k+1 = \ix k 
k Xk + A _1 w;fc(l + hu k+1 )(l + hv k+ i)x k+ i = fifa 



(5.39) 



Next, we transform the spectral equations further with the help of the evolutionary ones. 
From the first equations in ( |5.37|) , ( |5.39|) we derive: 



<t>k ~ (! + hv k )(l + hw k )<f> k = h\ii\) k 
and using this in the first equation in ( |5.39| ), we bring the latter into the form 

u k (l + hv k )(l + hw k )(f> k + hX{l + fru fc _i)(l + /iit» fc _i)0 fc _ 1 = /i(-0 fc - hX^ k _ x ) 
From the second equation in ( p.37[ ) we find the expression 



(5.40) 



1 + hw k )(l + hu k+1 )^ k+1 



1 



1 + hv k 

and plugging this into the second equation in ( |5.39|) , we find: 



(V'fc+i - hXip k ) 



1 + hv k 

Similarly, from the third equation in ( |5.37| ) we find the expression 

1 



(5.41) 



1 + hu k+1 )(l + hv k+1 )x k+1 



1 + hw k 

and plugging this into the third equation in (|5.39|) , we find: 



Xk+i ~ h ^Xk 



(5.42) 



Putting the equations ( |5.40|) , ( |5.41|) , ( |5.42| ) into the matrix notations, we finally find: 

r t~ 1 ua 2 a^> = $ 



{ A^W X = v4 
The compatibility conditions of ( |5.38| ), ( |5.43| ) coincide with ( |5.29|) . 



(5.43) 
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5.4 Explicit dBL2 -> " relativist ic" BL2 

Our general philosophy and the previous theorem suggest the matrix 

L + = F- 1 UA 2 WAz 1 V 



(5.44) 



as the Lax matrix of the "relativistic" BL2 hierarchy. The spirit of ||GK1|| , however, would 
suggest 

L + = T~ X UWV (5.45) 

and this is the choice we actually follow in the next Theorem. However, there is no conflict 
between these two suggestions, because changing 

W -> A 2 WA 2 X 



amounts solely to the simple change of variables 

1 + hv k 



w k -> w k 



1 + hv k+1 



(5.46) 



(not touching u k , v k ). In other words, in order to make the explicit discretization (|5.27| ) 
belonging to the "relativistic" BL2 hierarchy defined below, one needs to perform the change 
of variables ( |5.46| ). It turns out that the continuous time flows look more symmetric, if 
the Lax matrix is choosen as in (|5.45| ), and the discrete time equations ( |5.46| ) look more 
symmetric by the choice ( |5.44|) for the Lax matrix. 

Theorem 5.4 The hierarchy associated with the Lax matrix ( |5.45| ) is described by the fol- 
lowing Lax triads: 



U = AiU - UA 3 , V = A 2 V - VA , W = A 3 W - WA 2 

which have to be supplemented by the identity 

F = o = A X T - TAq 
so that for the matrix L + = J r ~ 1 UWV there holds the Lax equation 



Here 



L + — Aq, L + 
A = n4f(L + ))-a(f(L + )) 



(5.47) 

(5.48) 

(5.49) 
(5.50) 
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and the matrices -4.1,2,3 are uniquely defined by the condition of compatibility of the equa- 
tions of the Lax triads Q5.47| ), Q5.48| ). In particular, for the "first" flow of the hierarchy 
corresponding to f(L) = L, one has: 

N 

A) = ^ J2( w k-i v k + u k v k + u k w k + hUk^Wk^v^Ek^ 
fc=i 

N N 

+A~ 2 u k w k v k+1 E Kk+2 - hJ2 lkE kk (5.51) 
fc=i fc=i 

N 

Ai = \~ l Y( w k-i v k + u k v k + u k w k + hu k w k v k+ i)E k)k+1 

k=l 

N N 

+A~ 2 Y u k w k v k+1 E k>k+2 ~ hJ2 lk~\E kk (5.52) 

k=l k=l 

N 

M = \~ 1 Y( u k w k + v k w k + v k u k+1 + hu k v k w k )E ktk+1 

k=l 

N N 

+A~ 2 Y v k u k+1 w k+1 E Kk+2 ~hJ2 lkE kk (5.53) 

k=l k=l 

N 

A 3 = A -1 J2( v kU k +i + w k u k+1 + w k v k+1 + hw k u k+1 v k+1 )E k>k+1 

k=l 

N N 

+A~ 2 Y w k v k+ iu k+2 E k)k+2 - h Y lkE kk (5.54) 

k=l k=l 

where 

7 fc = w k ^v k + u k v k + u k w k + hu k _iw k _iv k + hu k w k v k+1 (5.55) 
The equations of motion of the "first" flow read: 

ii k = u k (v k + w k - - it> fc _i + hw k v k+1 - hw k - X v k + hj k - h^ k -i) 
v k = v k (w k + u k+1 - w k _i -u k + hu k w k - hu^iWk^ + hj k+1 - hj k ) (5.56) 
iv k = w k (u k+1 + v k+ i -u k -v k + hu k+1 v k+1 - hu k v k + hj k +l - hj k ) 

Proof is analogous to the proof of Theorem [O . We define Aq by (|5.50|) in accordance with 
the prescription of [|GK1|| , then put A\ = TA§T~ X G g ®g_. The matrices A 2 ,3 G g ®g_ are 
uniquely defined by the condition of consistency of the Lax triads in ( |5.47| ) . The results for 
the "first" flow ( |5.56|) (= "relativistic" BL2) are merely the specialization of this construction 
for f(L) = L and follow by a direct calculation. ■ 

The previous theorem allows to find the equations of motion of all flows of the hierarchy, 
and to calculate its integrals. Unfortunately, we still do not know its Hamiltonian structure, 
so that the problem of its integrability in the Liouville-Arnold sense remains open. 
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6 Conclusion 



We discussed in this paper an approach to finding new integrable hierarchies based on the 
time-discretization. Namely, performing the discretization of a known integrable lattice 
system and then changing the Cauchy problem for the discrete time equation, we arrive at 
another discrete time equation which belongs, generally speaking, to an integrable hierarchy 
disctinct from the one we started with. This new hierarchy is a regular deformation of the 
initial one, the deformation parameter being the time-step of the discretization. 

In the Toda lattice case the new hierarchy turns out to be the relativistic Toda one, the 
deformation parameter aquiring the meaning of the inverse speed of light. This suggests to 
call in general the new hierarchies obtained by this procedure "relativistic" ones. It would 
be interesting to find out, whether this notation allows some physical justification. 

In any case, we have now a program of finding "relativistic" deformations of a large 
number of known lattice hierarchies. In this connection the collection of local discretizations 
in [B6[ may serve as a starting point. For the systems belonging to the lattice KP hierarchy, 



another procedure of finding "relativistic" deformations was suggested in ||GK1|| , which has 



a priori nothing in common with our one. However, in the examples of the present paper, 
both procedures lead to similar results. It would be important to find out whether there are 
some deeper principles behind this coincidence. 
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